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H | Abstract. We study the dynamics of an asymmetric intruder in a glass- former 

model. At equilibrium, the intruder diffuses with average zero velocity. After an 
abrupt quench to T deeply under the mode-coupling temperature, a net average drift 
is observed, steady on a logarithmic time-scale. The phenomenon is well reproduced 
in an asymmetric version of the Sinai model. The subvelocity of the intruder grows 
with T e ff/T, where T e ff is defined by the response-correlation ratio, corresponding to 
a general behavior of thermal ratchets when in contact with two thermal reservoirs. 
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In an irreversible environment, thermal fluctuations can be rectified in order to 
produce a directed current. After a few fundamental examples of historical and 
conceptual value [II E], in the last twenty years a huge amount of devices and models 
— usually known as Brownian ratchets or motors — have been proposed [21 HJ [5] . The 
purpose of these models is often practical, e.g. the extraction of energy from a highly 
fluctuating environment, such as a living cell [61 |7]. But Brownian ratchets are also 
valid probes for the non-equilibrium properties of the fluctuating medium, the value of 
the current being sensitive to the interplay of different time-scales as well as different 
temperatures at work [9] . 

The simultaneous breaking of space and time-reversal symmetries is a necessary 
condition to observe directed motion. A classical example of this mechanism (see 
for instance [3]), the so-called flashing ratchet, consists of a particle undergoing 
standard overdamped diffusion and also subject to a zero average space-asymmetric 
(e.g. sawtooth) potential. If this potential is switched on and off according to a 
random/periodic time sequence, a current of particles (steady on average) can be 
observed. This current is driven by the energy injected into the system when switching 
on the potential. 

Another way of breaking the time-reversal symmetry is obtained by coupling the 
system with reservoirs at two different temperatures: an asymmetric intruder in such 
a multi-temperature environment displays an average drift, performing as a Brownian 
motor. In many-body systems this has been done, for instance, in [TUl [TT] . 

The present work, inspired from the latter scenarios with a continuous flow of 
energy in a non-thermalized medium, shows a study of the ratchet phenomenon in 
the aging dynamics of fragile glass-formers. It is well known that such systems, 
when quenched below their Mode- Coupling temperature, display an out-of-equilibrium 
dynamics customarily described within a two-temperature scenario [121 Ell O 115] - 
Fast modes are equilibrated at the bath temperature while slow modes remain at a 
higher effective temperature T e ff. Here we show that the energy flowing from slow 
to fast modes can be rectified to produce a directed motion. The properties of the 
observed current characterize the non-equilibrium behavior of the glass. In particular 
a striking monotonic relation is observed between the ratchet sub-velocity and T e ff/T. 
The experience with other kinetic models of ratchets [TTJJ [TT] teaches that - when in the 
presence of a temperature unbalance - the heat flux also governs the ratchet velocity. 
These observations suggest the conjecture that one can obtain the effective temperature 
by replacing the measure of linear response and correlations with the simple measure of 
an average current. In what follows we test the reliability of this procedure. 

The first of the numerical experiments proposed here involves the 3D soft-spheres 
model, which is a well known fragile glass-former [161 C3 CHI CHj • The thermodynamic 
properties of soft-spheres are controlled by a single parameter T = pT~ l l A , which 
combines the temperature T and the density p of the system, with p = N/Va^ and <t 
the radius of the effective one-component fluid. We study a binary mixture (50:50) with 
radii ratio, <Ji/a = 1.2. Particles interact via the soft potential U(r) = [(a.; + aj)/r} 12 



Ratchet effect in an aging glass 3 

and the dynamics is evolved via a local Monte Carlo algorithm. Time is measured in 
Monte Carlo steps (one steps corresponds to N attempted MC moves). The model has a 
dynamical crossover at a mode-coupling temperature Tmc corresponding to the effective 
coupling Tmc = 1-45 [16J. The ratchet is formed by the asymmetric interaction of a 
single particle, the intruder, with all the other particles of the system. Denoting with 
Xi, i > the x coordinate of the i-th particle, and with x the abscissa of the intruder, 
we choose 

JT( ,_) I7(|r — Ti|) ifx,<x , m 

1 °' l) \ eU(\r -ri\) otherwise, [ } 

with e = 0.02. The spatial symmetry along the x-axis is therefore broken, fulfilling one 
of the requirement to get a ratchet device. Let us stress that here, in analogy with 
the model proposed in [20], the asymmetry is inherent to a single object, the intruder, 
which is embedded in a symmetric environment, at variance with flashing and rocket 
ratchets [5]. Moreover, by exploiting condition (TfJ, we are able to build an intruder with 
an intrinsic asymmetry, even if it is a point-like particle. 

We equilibrate configurations embedding an asymmetric particle at a high 
temperature T iiq ~^> T MC - At T Kg the system has simple liquid behaviour with fast 
exponential relaxation. The dynamics of the asymmetric intruder is then studied both 
along equilibrium trajectories and after quenches to different temperatures T < T MC - 
The average displacements of the intruder along different axes is compared at and 
outside equilibrium, namely both (Ax (i)) and (Ay (t)) are studied, with Ax (t) = 
xo(t) —xo(to)- In particular t is the time elapsed since the instantaneous quench, namely 
we take to — 0. The averages denoted by (...) are realized considering 5000 initial 
configurations with a single intruder, equilibrated at Tu q = AA2T MC , each followed 
by an independent thermal history. Quenches are done to T = 0.67 Tmc, 0.53 Tmc, 
0.42 T MC , 0.31 T MC . 

Let us consider first the effect of the asymmetric interaction on the arrangement 
of particles around the intruder. We study the pair distribution function of neighbors 
to the right and to the left of the intruder, g^(r\Sx ^ 0) = Ylj\Sx >o^( r io — r ); with 
7jo the distance between the intruder and the j-th particle, in equilibrium or at a fixed 
elapsed time after the quench (see Fig. [H left panels). Right neighbors stay closer to 
the intruder due to the reduced repulsion. This asymmetric clustering of neighbors 
is slightly enhanced during aging, but at this level of analysis there is no qualitative 
differences between the equilibrium and off-equilibrium regimes. 

Dynamical measurements, on the contrary, show important differences between 
the two regimes. Let us consider the mean square displacement around the average 
position at time t of the intruder particle (Fig. [1] right): it is linear with time at high 
temperatures, while it is logarithmic after a quench, (Axq(£)) — {Ax (t)) 2 ~ logt. 
This behavior, typical of activated dynamics in a rough potential J2TJ [T71 [22], is 
also observed for host particles. The behavior of the average displacement of the 
intruder, (Axo(t)), is more striking. At equilibrium (Fig. [21 black circles) there is no 
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Figure 1. Left panels: pair distribution function centered on the intruder, g^(r\5x ^ 
0), of left (red squares), and right (black circles) neighbors, at equilibrium (top) and at 
a fixed elapsed time after quench (bottom). Black line, g(r) for symmetric interactions. 
Right panels: (A:tq(£)) at equilibrium (top) and after the quench 

(bottom). 



net displacement: this is because parity along x axes is broken whereas time reversal 
symmetry is preserved. The same happens to (Ay (t)) (green diamonds) during aging: 
in this case only (macroscopic) time reversal symmetry is broken. But parity and time- 
reversal symmetry are both violated for (Axo(t)) after the quench (red squares), and 
in this case a net average drift is found (Fig. [2]), linear on a logarithmic timescale, 
(Axo(t)) ~ r with r = log ' t. The logarithmic timescale again points to a non- 
equilibrium phenomenon ruled by activated events. Moreover, the simple scaling relation 
(Axo(i)) ~ a/(Axo(£)) — (Axo(t)) 2 between the displacement along the x axis and the 
m.s.d. around the average position is observed, in analogy with the Sinai model discussed 
below. 

Our first exploration of ratcheting effects in glassy models leaves open one important 
question, namely the dependence of the drift on the parameters (e.g. temperature) of 
the system. In order to sketch a preliminary answer, we need to extract from the curve 
(Axo(t)) a synthetic observable. The finding of a logarithmic timescale (Axo(t)) 
with r ~ log ' t suggests to define an average sub-velocity as [23] : 

(Ax (t)) - (Ax {t w )} (x (t)) - (x (t w )) 



T 



Vsubyt, t t 



(2) 



5t 5t 

with 5t = log x ' 2 t — log 1 ' 2 ^. This average subvelocity depends in general on both 
the running time t and the waiting time t w elapsed since the quench. More precisely, 
considering fig. [2j by fixing t w and t we choose the time lag where the slope of the 
curve (Axo(t)) is measured. Clearly, the instantaneous sub- velocity only depends on t w 
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Figure 2. Average intruder displacement: (Acco(i)) for equilibrium trajectories (black 
circles), and (Axo(i)) and (Ay (t)) after a quench far below Ta/c {T = 0A2Tmc) (red 
squares and green diamonds respectively). The displacement (Axo(t)) is measured in 
units of the average inter-particle distance, as obtained from the position of the first 
peak in the pair distribution function g(r) (see also fig. [J). 



and corresponds to a slowly decaying velocity v(t w ) = d(Ax (t w ))/dt w ~ l/t w (with 
logarithmic corrections). For large enough waiting times we can define the "order 
parameter" v su b, namely we find v su b(t w ) ~ const, and then we probe its dependence 
upon the external parameters. Here we focus on the quench temperature T, drawing 
a connection between the drift of the glassy ratchet, which is a pure non-equilibrium 
effect, and the more customary equilibrium-like descriptions of aging media in terms of 
effective temperatures. In the right panel of Fig. [3] we show the drift of the asymmetric 
intruder for four quench temperatures: T = 0.67T MC , 0.53T MC , 0A2T MC , 0.31T MC - 
We observe that when decreasing the quench temperature the drift grows in intensity, 
i.e. the sub-velocity increases. This somehow is a stronger evidence that the ratchet 
drift cannot be described as an equilibrium-like effect, e.g. trying to connect average 
kinetic or potential energy with a mobility: one would expect that at lower temperatures 
everything is slowed down, whereas we find that the ratchet drift is enhanced. 

Typical examples of ratchets in (ideally) statistically stationary configurations are 
obtained by coupling the system with two or more reservoirs. It is the existence of 
different temperatures within the same system which allows the production of work 
without violations of the second principle of thermodynamics. But what is the second 
temperature in a glassy system? According to the well-established description of the 
aging regime of glasses p3], it is the effective temperature defined as the violation 
factor of the fluctuation-dissipation theorem (FDT). The quench of a fragile glass, for 
instance our the soft spheres model, below its mode-coupling temperature produces 
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aging and violations of FDT. The factor X(t,t w ) < 1 which allows one to write a 
generalized FDT, Tx(t,t w ) = X(t,t w )[C(t,t) — C(t,t w )], with x(^>0 the integrated 
response and C(t,t w ) the correlation, yields the definition of the effective temperature 
T e ff(t,t w ) = T/X(t,t w ). The last is usually higher than the bath temperature 
T e ff(t,t w ) > T and is understood as the temperature of slow, still not equilibrated, 
modes. Clearly, the ratio T e ff/T may be regarded as the parameter which tunes non 
equilibrium effects and we study here how the glassy ratchet drift depends on it. We 
obtain T e ff from the parametric plot of C(t, t w ), taken as the self-intermediate scattering 
function, versus the integrated response Tx(t,t w ), measured according to the field-free 
method of [21] . The parametric plot Tx{t,t w ) vs C(t,t w ) is shown in the left panel of 
fig. [3] for different temperatures. 
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Figure 3. Left panel: parametric plot Tx(t,t w ) vs C(t,t w ). Right panel: asymmetric 
intruder drift for different quench temperatures; linear fits yield sublinear velocities. 
Inset: v 3U t, vs T e ff/T. In both panels colors of data correspond to different effective 
tempeartures: T eff = 1.497 T MC (black), 1.613 T M c (red), 1.792 T M C (green) and 
1.819 T M c (blue). 



At large times v su b exhibits finite-size effects, namely the drift saturates at a time 
that can be increased by increasing the size of the simulation box. This is why we 
compute v su b by fitting data only in a relative early time window, namely 10 < t w < 10 3 . 
Accordingly, the measure of X(t,t w ) is obtained from the parametric plot Tx(t,t w ) vs 
C(t,t w ) with t = 10 3 and t w ranging from 10 to 10 3 , namely the early aging regime is 
considered. The inset of the right panel of fig. [3] shows the behavior of v su b vs T e ff/T, 
revealing that the subvelocity increases when T e ff/T is increased. Namely the intensity 
of the ratchet effect, traced in the measure oiv su f,, grows as the distance from equilibrium 
is increased. 

Let us remark here that the mechanism governing our ratchet differs from that of a 
flashing ratchet [3]. Irreversibility is achieved by choosing an initial condition which is 
out-of-equilibrium with respect to the bath temperature Tu q ^ T. The extreme slowness 
of an aging glass prevents the system from thermalizing, so that energy continuously 
flows from fast to slow modes, supplying power to the Brownian ratchet. At variance 
with a flashing ratchet, here there is no external time- dependent modulation of a 



Ratchet effect in an aging glass 7 

potential. The movement of surrounding molecules produces fluctuations of the intruder 
potential, but these obey a globally conservative dynamics and cannot produce energy. 

Indeed, simpler models can be conceived, where the intruder moves in a time- 
independent potential: as a enlightening example, in the following we discuss the Sinai 
model, where it is made clear how irreversibility comes solely from the choice of initial 
conditions. Sinai model is one of the simplest describing the diffusion of a single particle 
through a random correlated potential [25]. Its long-time dynamics is ruled by activated 
events and is characterized by a logarithmic time-scale: for this reasons it appears to 
be a well fitted candidate to reproduce the previous experiment in a more controlled 
setup. In the original Sinai model the random potential is built from a random-walk of 
the force on a Id lattice. The force F^ at each lattice site i is an independent identically 
distributed random variable extracted from a zero mean symmetric distribution p(F). 
The potential on a lattice site n is given by U(n) = Y^i=i F%- Because f dFp(F)F = 
the average force experienced by a particle is zero. The potential excursion between 
two sites grows like (\U(i) — U(j)\) ~ \i — j] 1 ^ 2 . The above relation, together with the 
expression of the typical time needed to jump a barrier, ~ exp(/3A£7), yields the long- 
time scaling of the mean square displacement of: (Al 2 (t)) ~ log (t), with l(t) denoting 
the site occupied by the particle at time t and Al(t) = l(t) — 1(0). A logarithmic time 
scale for the growth of domain size is quite ubiquitous in the low temperature regime of 
glassy systems, where activated processes dominate [26J. The time- reversal symmetry 
breaking is inherent in the Sinai model: by averaging over initial positions extracted 
from a flat distribution one reproduces an initial infinite temperature Tu q = oo, while the 
quench to a glassy phase is reproduced evolving the system at temperature T <C y/L 
where L is the linear size of the system. We propose a spatially asymmetric version 
of the Sinai model, where the symmetric force distribution p(F) is replaced with an 
asymmetric one p(F), in analogy with the asymmetric potential of the glassy ratchet, 
in order to find a non-equilibrium drift. Namely the distribution p(F) is such that 
p(F) ^ p(—F), but still J dFp(F)F = 0. In particular, the random force is obtained 
according to the following procedure: at each site the sign of a random variable fi is 
chosen with probability 1/2 and its modulus, in agreement with the sign, from either 
e -|/*IM- or e -|/»l/A+ j with A + > A_. In the present simulation A+ = 1 and A_ = 0.2 have 
been used. The potential is then built as U(n) = Y^i=\ -^> with Fi = f — (A+ — A_)/2. 
This amounts to a shift of the whole distribution such that (F) = 0. 

The results of the simulations of this model are shown in Fig. HI they clearly show 
a behavior in striking similarity with those of Fig. [2] for the glass-former model: a drift 
is observed only when both time reversal and spatial symmetry are broken. Indeed, 
as demonstrated by the black and red solid curves, it is sufficient that one of the two 
symmetries is restored to have a zero drift. In particular, to obtain an "equilibrium" 
dynamics in this model it is sufficient to distribute all initial positions of the simulation 
according to a distribution ~ exp[—/3U(i)]: in this case, even with the asymmetric 
distribution of forces described above, the average drift is zero. The large time behavior 
of the drift is compatible with a squared logarithm, (Al(t)) ~ log 2 (t) ~ ^/(A/ 2 (t)), in 
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Figure 4. Main: Average displacement in Sinai model. Black circles, diffusion 
in asymmetric potential at equilibrium; green diamonds, diffusion during aging in 
symmetric potential (quench from T = oo); red squares, diffusion during aging in 
asymmetric potential. Inset: the sublinear drift velocity v su b grows with \JLjT . 



agreement with what observed previously for the glassy ratchet. 

The out-of-equilibrium dynamics of the Sinai model closely reproduces these 
observations: in the inset of Fig. @] shows the sub- velocity as a function of a/Z/T, with 
T <C L 1 / 2 the quench temperature and L the size of the linear chain, which fixes the 
most relevant energy scale of the model, i.e. the maximum depth of potential minima. 
The inset of Fig. H] shows that v su b grows monotonously with yLjT: when the last 
quantity is increased particles condensate on the bottom of the deepest valleys. The 
monotonic increase of v su b with vL/T signals that also for the Sinai model the larger 
the distance from equilibrium the larger the velocity of the drift. 

In conclusion, through numerical simulations in different models and different 
choices of the quench temperature, always chosen in the deep slowly relaxing regime, 
we have given evidences of the existence of a "glassy ratchet" phenomenon. The drift 
velocity slowly decays in time and can be appreciably different from zero for at least 
three orders of magnitude in time. The overall intensity of the drift, measured in terms 
of a "sub- velocity" , is monotonically increasing with the distance from equilibrium, i.e. 
with the difference between the quench and effective temperatures. This observation 
supports the idea of regarding the ratchet drift as a "non-equilibrium thermometer" : it 
can be used as a device capable to say how far is a system from equilibrium. Nevertheless, 
for such a thermometer to be effective, a more accurate calibration procedure should be 
carried on. Namely, one should verify that v su b is a function of T e ff and T only with 
a small number of parameters. In this case the calibration of our thermometer would 
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require a small independent measures of T e ff to fix those parameters. 

As an experimental realization of our glassy ratchet, one should consider a particle 
with anisotropic interaction with the surrounding ones with a small magnetic dipole 
placed on it, orthogonal to the asymmetry axis. Then, switching on a constant magnetic 
field the orientation of the two faces of the "Janus" [27] particle is preserved with respect 
to a fixed reference frame, so that the spatial symmetry of the interaction is broken. 
Recent theoretical and experimental advances in the study of functionalized or "patchy" 
particles [28j EH] promise an experimental verification of our hypothesis in the near 
future. 
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